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Polymer brushes are increasingly used to tailor surface physicochemistry for various applications 
such as wetting, adhesion of biological objects, implantable devices, etc. We perform Dissipative 
Particle Dynamics simulations to study the behavior of dense polymer brushes under flow in a slit- 
pore channel. We discover that the system displays flow inversion at the brush interface for several 
disconnected ranges of the imposed flow. We associate such phenomenon to collective polymer 
dynamics: a wave propagating on the brush surface. The relation between the wavelength, the 
amplitude and the propagation speed of the flow-generated wave is consistent with the solution of 
the Stokes equations when an imposed traveling wave is assumed as boundary condition (the famous 
Taylor’s swimmer). 

PACS numbers: 47.56.+r, 47.63.-b, 82.35.Lr, 87.15.H- 


Polymer brushes are passive media whose great variety 
allows for a rich range of applications. Brushes with dif¬ 
ferent mechanical properties can be created by grafting 
simple polymers, block copolymers or polymer stars to 
a solid substrates or to an interface between two liquids. 
They are exploited for colloid stabilization m, as lu¬ 
bricant layers BH2], as adhesion regulators EH] and for 
biomedical technological applications. A holdover inter¬ 
est in these systems is indeed motivated by the discovery 
that the inner surface of various mammalian organs is 
decorated by densily grafted macromolecules. For exam¬ 
ple, the lumen of blood vessels is protected by a hun¬ 
dred nanometers thick polymer brush mainly made of 
glucose and its compounds. Such an heterogeneous net¬ 
work is called “glycocalyx” mm- Research on the gly¬ 
cocalyx dynamics is central for a complete understanding 
of the blood circulatory system and, hopefully, to shed 
light on the initial stages of diseases like thrombosis and 
atherosclerosis m- Understanding glycocalyx is also 
relevant to the polymeric coatings of lungs, small intes¬ 
tine and uterus H2HH]- 

Most theoretical studies have attempted to model 
polymer brushes as porous media described by 
Brinkman-type equations m or as elastic media com¬ 
posed of very rigid fibers M- In this Letter, we report 
an analysis of a polymer brush whose features approach 
the glycocalyx system and for which the set up conditions 
recall the microcirculation frame na. We focus on dense 
polymer brushes under flow, in which the brush is mod¬ 
eled as a collection of individual polymers. We implement 
a Dissipative Particle Dynamics (DPD) code [TM2TI with 
explicit solvent and numerically analyse the dynamics of 
a linear flexible honro-disperse polymer brush subdued 
to a simple liquid parabolic flow in a slit-pore geometry. 
The coarse-grained DPD procedure applies to both sol¬ 
vent molecules and polymer monomers, offering the pos¬ 


sibility to (i) reproduce hydrodynamic interactions while 
retaining a detailed view of the brush dynamics on the 
scale of the coarse-grained monomers; (ii) access both the 
polymer dynamics, influenced by the imposed flow, and 
the flow field, perturbed by the presence and motion of 
the brush. 

Recent studies of polymer brushes under flow [?3 - lMl 
have highlighted an unexpected behavior in the veloc¬ 
ity profile in the vicinity of the brush surface. These 
studies have reported that the velocity field reverses on 
increasing the flow field and have tentatively associated 
such result to the peculiar dynamics of the single poly¬ 
mer undergoing a cyclic motion of stretching, elongation 
and recoiling. Stimulated by these studies we have un¬ 
dertaken a numerical investigation exploring a very large 
range of imposed flow, aiming at quantifying the condi¬ 
tions under which the physical properties of the brush 
couple with the hydrodynamic properties of the solvent 
to produce flow inversion. As presented below, we dis¬ 
cover (i) that flow inversion appears in distinct regions of 
imposed flow values; (ii) that every time flow inversion 
is observed a surface wave appears, stressing that such 
backflow is strongly associated to a collective (as opposed 
to single) polymer dynamics; (iii) that the wave proper¬ 
ties are consistent with predictions derived by Taylor in 
his seminal study on pusher microswimmers |25j . Thus, 
our work presents a new interpretation for the flow in¬ 
version phenomenon and provide a novel connection be¬ 
tween two separated fields: polymer brushes under flow 
and microswimmers. 

A short description of the DPD methodology and of 
its application to the slit-pore geometry is reported in 
the Supplementary Informations (S.I., Sect.l). Here we 
recall that we simulate a parallelepiped box of sides L x , 
L y , L z . with polymers chains composed by 40 (coarse¬ 
grained) monomers grafted at z = 0 (see Fig. [I]). Our 
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FIG. 1. Schematic representation of the investigated slit-pore 
system. Polymers are grafted at the bottom wall z = 0 while 
the rest of the channel is occupied by solvent particles. An 
additional sketch is provided with a zoom on a typical chain 
(highlighted in red). The constant acceleration A is applied 
along the x axis to all solvent particles. 


discussion pertains simulations in a box of size L x = 30, 
L y = 5 and Lz = 50, in units of r c (see S. I. for defini¬ 
tions). We remark that we observe the same phenomenol¬ 
ogy we are going to present also in considerably bigger 
boxes ( L x = 360, L y = 20 and L z = 150). However, for 
the sake of computational time and memory capacity, 
the systematic analysis has been conducted in a smaller 
channel. The equilibrium distance between neighbour 
monomers in a chain is 0.89r c . The grafting density 
cr graft = N ch /(L X L V ) with N c h the number of chains 
composing the brush, being cr gra f t = 1.5, corresponds to 
the condition of dense brush. Periodic boundary con¬ 
ditions are applied along the x and y directions. In the 
channel, a flow is imposed by applying a constant acceler¬ 
ation Ax to all solvent particles, resulting in a parabolic 
velocity profile along the z direction. The strength of the 
velocity is controlled by the value of A. In the follow¬ 
ing, instead of A , we will use the so-called Weissenberg 

number Wi = brush ; which provides an adimensional 

^ flow 

information on the relative timescale of the unperturbed 
brush dynamics ( tb rus h ) over the inverse of the averaged 
shear rate inside the channel (tfi ow ). Specific values for 
all parameters and their units and a detailed discussion 
of the A dependence of tfi ow are reported in the S.I., 
Sect.1,3. 

The velocity profile v x (z), calculated by averaging the 
^-component of the velocity for each solvent particle as 
a function of its position z, is shown in Figs. H(a)-(b) 
for two Wi values, providing two typical examples. As 
expected, the functional shape of v x (z) in the region 
around the largest velocity v ma x is well represented by 
a parabolic function with the expected amplitude. In all 
cases, the parabolic function predicts that v x vanishes 
when 2 approaches the brush height h, suggesting that, 
for all studied Wi, the presence of a dense polymer brush 
restricts the pore by an amount equivalent to h. The ve- 



FIG. 2. Panels (a) and (b) show the velocity profile in the 
region where the brush is set and they are a zoom of the 
respective insets. In both cases, corresponding to (a) Wi = 45 
and (b) Wi = 418, the profile is parabolic around Vmax■ In 
(a) no backflow is registered, while in (b) the flow inversion at 
z ~ h is evident. Panel (c) reports the minimum value v m in 
of the v x (z) profile for different Wi (symbols). The line is a 
guide for the eye. It is possible to recognize three regions for 
which Vmin assumes a negative value, signalling the presence 
of a backward flow. 


locity profile in the region z < h, i.e. the region occupied 
by the brush, is particularly interesting. In some cases, 
e.g. the one of Fig. [2] (a), the fluid velocity inside the 
brush is small, consistent with a picture in which the hy¬ 
drodynamic interactions are effectively screened by the 
presence of the polymer layer. In other cases, as shown 
in Fig [2j(b), the velocity profile exhibits a flow inver¬ 
sion at the interface between the brush and the bulk, e.g. 
around z ~ h, as previously documented for an analo¬ 
gous system JJj as well as when the solvent is replaced 
by a polymer melt in [22[G3]- In these previous studies, 
the onset of flow inversion was associated to the point at 
which the shear rate exceeds a certain threshold )24j . Op¬ 
positely, we find that flow inversion takes place in differ¬ 
ent distinct disconnected windows of Wi values. Indeed, 
we uncover at least three different windows of flow inten¬ 
sities in which a backflow is observed. In between these 
regions the velocity profile at z ~ h resembles the one in 
Fig.§(a). Fig.§(c) shows the minimum value v m i n of 
the profile v x (z) near the surface of the polymer brush, 
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as function of Wi. By this definition, flow inversion is 
observed when v m i n < 0. The identified three back-flow 
regions fall around Wi ~ 200,400, 600. 

In the previously mentioned papers it was ar¬ 

gued that the backflow is a consequence of the grafted 
single chain dynamics, performing a recursive imperfect 
cycle composed of: tilting, elongation and recoiling. In¬ 
deed, an elongated chain (and thus exposed to flow) is 
dragged by the shear stress along the flow direction and 
then recoiled back by entropic forces. By examining the 
polymer trajectories we confirm the presence of such re¬ 
cursive motion (see movie in the S.I., Sect.4) for all Wi 
values, but such cyclic motion is found to be indepen¬ 
dent on the presence of the back-flow effect. Hence, such 
single-polymer motion can not explain the onset of flow 
inversion. 



(a) (b) 



(c) 


FIG. 3. Sketches of the whole brush, corresponding to (a) 
Wi = 91 and (b) Wi = 418. In (a) the brush surface is basi¬ 
cally flat, while in (b) the surface is modulated by a traveling 
wave with \ ~ L x = 30, v = 0.077 and b = 1.18. In (c) a 
snapshot for a larger studied systems (L x = 180, L y — 5 and 
L z = 150) with a surface wave with multiple crests of wave¬ 
length L x /2 « 90 (Wi = 340). We show the central box (in 
red) and the two nearest replicas (in grey) to highlight the 
periodic nature of the waves. 


Investigating space and time correlations among dif¬ 
ferent polymers, we discover a wave traveling over the 
brush surface in the same direction as the imposed flow. 
The wave arises only when back-flow is present, suggest¬ 
ing a strong association between the surface wave and 
the capability of the brush to produce an inversion of 
the flow velocity around z ~ h. A visual observation 
of the brush and of its dynamics nicely show the col¬ 
lective behavior of the polymers: Fig. |3](a),(b) show a 
typical frame for two generic cases of absence and pres¬ 
ence of flow inversion. Videos of the time dependence 
of the brush are available in the S.I., Sect.4. We find 
that this traveling wave has a clear spatial periodicity 
of the order of the simulation box and a non-negligible 
amplitude. We characterize such surface wave evaluating 
its frequency u, wavelength A and oscillation amplitude 
b. To do so we define the brush surface S(x,t ) as the 
position of the furthest monomer from the grafting wall 
occupying at time t the coordinate x, averaged over all 
y. The power-spectrum of the time Fourier transform of 


S(xo,t ) (for an arbitrary xo value, inset of Fig. @(a) ) dis¬ 
plays, in all cases in which a flow inversion is observed, 
a sharp peak at a given frequency (Fig. @( a)) that we 
associate to v. In some cases a much less intense peak 
is also observed at 2za Fourier transforming the signal 
S(x,to) in real space for an arbitrary to and averaging 
over all configurations provides a quantification of the 
wavelength A. As shown in Fig. |4](b), for all cases in 
which a flow inversion is present we find, as in the time 
domain, a dominant contribution from the longest A that 
can propagate in a period system of size L x and from its 
second harmonic (e.g. from wavenumbers k equal to 
and j^-). Fig. Ftl(c) shows the Wi dependence of v for 
systems with the same L x (e.g. same A). In each of the 
three disconnected regions of Wi values where flow inver¬ 
sion is present, v increases approximative^ linearly with 
Wi. Finally, we define the wave amplitude b as the stan¬ 
dard deviation of S(x,t) and the wave average height as 
S=< S(x,t) >. 

We stress that wave propagation is observed also in sys¬ 
tems of different L x . We find that increasing the length 
L x and the width L z of the channel, multiple wave crests 
can be displaied (in Fig. i( c) an example). This sug¬ 
gests that the finite size of the simulation box exerts a 
cutting on the density of states, making it possible to 
observe only waves with wavenumbers which are integer 
multiple of 2-k/L x . Further details on finite size effects 
are reported in the S.I., Sect.5. 

To dig into the reasons why a flow inversion is ob¬ 
served every time there is a wave propagation, we recall 
the derivation provided by Taylor [25], dated back to 
1951, to explain the mobility of pusher microswimmers. 
In his seminal paper, Taylor determines analytically the 
self-propelling velocity Vrayior, at low Reynolds number, 
of an infinite two dimensional sheet whose profile is a 
sinusoidal wave traveling with (opposite) velocity U: 


\V T ay 


2 7T 2 b 2 f 19 TT 2 b 2 \ 

A 2 \ 4 A 2 ) ’ 


(1) 


where b and A are, respectively, the amplitude and wave¬ 
length of the sinusoidal wave. By associating the two di¬ 
mensional sheet with the brush surface, the wave velocity 
of the sheet with the wave velocity of the brush surface 
and the self-propelling velocity of the sheet with the op¬ 
posite velocity of the solvent at the brush surface u m i n , 
Taylor’s expression provides a precise relation between 
the parameters entering in the wave phenomenon and 
the parameters controlling the flow inversion. In Fig. [5] 
we show, for each set of b , A and U values associated to 
a specific Wi value giving flow inversion, both \Vr a yior\ 
(from Eq. 0) and \vmi n \. All the trends of the \v m i n \ 
behavior are perfectly recovered by the self-propelling ve¬ 
locity of a Taylor’s micro-swimmer, suggesting that in¬ 
deed the elastic brush can be considered as an anchored 
micro-swimmer, unable to move but able to propel the 
surrounding liquid. 
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FIG. 4. Fourier power-spectra: (a) in the frequency domain 
u for the signal S(xo,t), (b) in the wavenumber domain k = 
2n/X for the signal S(x, to) (red squares indicate the spectrum 
in a bigger system L x = 180, L y — 5 and L z = 150). In 
both cases the first and second harmonics are visible. The 
two insets show typical signals in real space. Panel (c) shows 
the first harmonic frequencies corresponding to the Wi giving 
backflow. Such frequencies monotonically increase. 


As a last consideration, we investigate possible origins 
of the surface wave. Our results suggest that the wave 
propagation arises from hydrodynamics interactions syn¬ 
chronising polymers, in a sort of metachronal motion. A 
“metachronal wave” is known to develop in ciliated sys¬ 
tems, namely flat surfaces covered by a regular matrix 
of equally spaced flexible filaments pa 123- Exposed to 
flow, filaments do not interact directly, but thanks to hy¬ 
drodynamics interactions, after some transient time, they 
reach a synchronized state with a regular phase shift be¬ 
tween subsequent filament rows. However, in those cases, 
the matrix is composed by active matter elements and 
the sequence of configurations they assume is assigned a 
priori. On the contrary, polymer brushes are passive me¬ 
dia. Different approaches have been developed to model 
wave formation on passive media such as crop canopies 
(because of flowing wind) or aquatic vegetations (because 
of water flow), but they require also the inertial term of 
Navier-Stokes equations to account for the instability ini¬ 
tiating the surface modulation [251125] . However, we can 



FIG. 5. The mapping between the absolute value of the ve¬ 
locity profile minimum \v m in\ (stars) and the velocity Vrayior 
obtained by inserting in Eq. 0 the wave amplitude b, the 
wavelength A and the speed velocity U = Xu evaluated for 
the brush surface wave at the given Wi. 

consider the brush to be an elastic medium. Therefore, 
the presence of a wave suggests that a possible resonance 
effect may arise between the shear rate 7 at the brush sur¬ 
face, produced by the imposed flow, and the frequency 
u com associated to the time scale of relaxation from the 
elastic compression of the whole brush 130]. This would 
point to one threshold value of Wi beyond which surface 
oscillations are observed, which is not consistent with the 
finding of waves in at least three different Wi regions. To 
clarify this inconsistency, we show in Fig.[ 6 ]the shear rate 
7 estimated by a linear fit of the velocity profile v x (z) at 
the brush surface 2 : ~ S, as a function of A. We ascribe 
the striking non-linear behaviour of 7 vs. A to the com¬ 
pression of the brush at high imposed flow. We also find 
a significant dependence of the brush elastic compression 
properties on A, suggesting the possibility of multiple in¬ 
tersections of the two quantities 7 and u com (and hence 
multiple resonances and multiple regions of flow inver¬ 
sion) on varying A or, equivalently, Wi (see S.I., Sect.2). 

In summary, we have discovered for dense polymer 
brushes under flow a well-characterized surface travel¬ 
ing wave, which appears in all cases of back-flow and it is 
not present in the other ones, providing a novel picture 
of the flow-inversion phenomenon, associated to a collec¬ 
tive motion of the whole brush rather than to a single 
chain dynamics. We have also observed a striking and 
unexpected similarity in the relation between the wave 
velocity and the back-flow velocity in the case of flow in 
passive polymer brushes and the same relation for the 
case of active pusher microswimmers [25]. 

Discussions with prof. Philippe Peyla are kindly ac¬ 
knowledged. C.M. and S.B. also acknowledge finan¬ 
cial support from ONES (Centre National d’Etudes Spa- 
tiales) and ESA (European Space Agency). L.R. ac¬ 
knowledges support from the Austrian Research Fund 
(FWF) through his Lise-Meitner Fellowship M 1650-N27. 














5 



FIG. 6 . The shear rate 7 measured from a linear fit of the 
velocity profile v x (z) at the brush surface z ~ S for different 
flow intensities A. 
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Supplementary Information 


1. DPD METHODOLOGY 

DPD is a coarse-grained Molecular Dynamics method. 
It was introduced by Hoogerbrugge and Koelman in 
1992 [TS] to simulate isothermal Navier-Stokes equations 
and it has been recently recognized that it can be success¬ 
fully applied to any mesoscopic scale [31]. In the spirit 
of DPD, the simulated particle does not correspond to a 
single molecule, but rather to a significant large cluster 
of them. Each of the N point-like particles evolves in 
time according to Newton’s equations 

N N 

mj i = Y,F ij =J2( P !j+ P Z+ P & ( 2 ) 

that we solve by using the velocity Verlet algorithm [32] ■ 
Observables are then calculated as averaged all over time 
configurations. The force F,j of Eq. [ 2 ] on each particle 
has three contributions m resulting from the coarse- 
graining procedure: a conservative component Ffj. a dis¬ 
sipative component FR and a stochastic one F$. All 
these forces are pair-wise, to guarantee momentum con¬ 
servation, and limited to a cut-off radius r c . 

The conservative soft-core repulsive force has the fol¬ 
lowing expression: 

pc _ f a «>/3 (l — rij < r c , (gN 

\o rij > r c 

with fij = fi — fj is the vector distance between the i- 
th and j-th particle, ry,- = \ fij\ and r^ = fij /rij. Such 
repulsive potential models the averaged fast microscopic 
length and time scales and was first derived by coarse- 
graining particles interacting via a Lennard-Jones poten¬ 
tial [33]. The indexes a,/3 indicates the particle type 
(solvent or polymer, in our case). The constant a a> p 
measures the force between two completely overlapping 
particles and it is casted from the compressibility of the 
modeled fluid. 

The two other forces account for the loss of details in 
the coarse-graining procedure. The dissipative force has 
the form: 


p ij = -l w ° ( Uj)(hj ■ Vij)fij , (4) 


where the standard choice in literature for the “weight 
function” w D (r^) is: 


w D (nj) 



r^ < r c , 
Uj > r c . 


(5) 


Equation Q introduces a friction among particles pro¬ 
portional to the relative velocity = Vi — Vj Thermal 
fluctuations are added via the stochastic force: 


Ffj =aw s {rij)9ij(At) ( 6 ) 


where a is a constant (related to temperature), w s (rij) 
is a weight function and Otj is a random number ex¬ 
tracted from a gaussian distribution with zero aver¬ 
age and indipendent in time and among particle pairs: 
(' Oij(t)0im{t')) = {Su6jm + SimSji)S(t - t'). Momentum 
conservation requires 9ij = Oji. 

Two additional constraints, w D {rij) = [w R (rij)] 2 and 


7 = 


2 k B T 


, guarantee that the probability distribution 


respects the statistics of the NVT ensamble [20] . 

To mimic the brush chains we add a finite extensible 
nonlinear elastic potential (FENE) [3H 


F/f ne = -2 kR 2 


R 2 ~ ( r b _ r eq)' 2 


[j r e q < R (7) 


which account for the neighbour particle connectivity. In 
Eq. 0 , r eq is the equilibrium distance between neighbour 
monomers and R is the maximum allowed extensibility. 
Grafted points are randomly chosen from a uniform dis¬ 
tribution and located on a flat surface. Note that poly¬ 
mers are non-ideal since they interact via an excluded 
volume potential. 


1.1 DPD and physical units 

We chose the DPD units such as r c = 1, mi = 1, 
Idpd = 1 and ksT = 1. Following HD, we fix the 
solvent-solvent interaction parameter ass relating it to 
the adimensional compressibility of water ac^ 1 . A com¬ 
parison between the equation of state of a simulated DPD 
system and the experimental data (n^ 1 (water) = 16 at 
room temperature) suggests ass = 75 ksT/p = 25. We 
assume that the polymer-polymer interaction parameter 
has the same value (e.g. ass = app), while we se¬ 
lect a smaller value for the solvent-polymer parameter 
asp = 20 (good solvent conditions). The noise ampli¬ 
tude is fixed to er = 3. For the FENE potential, finally, 
we use r eq = 0.86, R = 1 and k = 50. 

In DPD the choice for the physical units depends 
on the level of coarse-graining is desired. Our refer¬ 
ence system is the endothelial glycocalyx, therefore we 
fix the physical lengthscale l p h ys relating the spacing 
between different filaments of the glycocalyx network, 
dgiyco = 20 nm [35; with the average distance between 
anchor points of our brush, d gra ft = ^/1 /ct graft = 0 . 82 : 

lphys dgly C Q/dg ra ft 24 • 10 TYl. 

We underline that with this choice also the brush thick¬ 
ness is in the range of endothelial glycocalyx (lOOnm 
-fTOOOnm [55]). 
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As for the physical mass scale m p h ys and physical time 
scale tphysi we exploit the comparison between viscosi¬ 
ties and between energies. The viscosity of water r/phys 
is Vphys = 10 -3 Pa • s at 300AT and the DPD viscosity 
Pdpd in a bare channel, estimated from the slit-pore 
velocity profile relation 77 = pA(L z /2) 2 /2v max , where 
Vmax is the maximum velocity (see Sect. 1.2 for defini¬ 
tions of the other parameters), corresponds in our case 
to r] DPD = 0.84. The physical energy scale is defined as 
ksTphys, with Tphys = 300AT. Thus we can write down 
two relations 


"Iphys 


r lphys ^ , 

l 'phys lj phys ■> 

VDPD 


t 


phys — 


Lphys 

Vphys 


Lphys 


3 kg Tphys 
YTlphys 


from which we extract 


tphys — 


/ 3 

l phys 


Vphys 


SksTphys 1]DPD 


= 1.8 • 10 _ 6 s, 


TP phys — 


Pphys + , 

Lphys^phys 

VDPD 


5.1 • 10 " 17 I<g. 


1.2 The system geometry 

We simulate a parallelepiped box of sides L x = 30, 
L y = 5 and L z = 50. On the x- and y- axes we im¬ 
pose periodic boundary conditions, while at z = 0 and 
z = L z two impenetrable parallel walls of infinite mass 
are set. At the wall we use the so-called “bounce back 
reflection” conditions m in which all the velocity com¬ 
ponents are reversed: 


v —> — v at z = 0,L z . (10) 

Such conditions assure no-slip boundaries. The integra¬ 
tion time step is At = 0.02. We set the density p = 
N/(L x L y L z ) = 3 |T5], therefore N = 22500. We fix the 
number of monomers per chain n = 40 and the grafting 
density cr gra ft = 1-5, defined as cr gra ft = N ch /(L x L y ), 
with N ch the number of chains composing the brush. For 
the sake of ease and computational time, we attach poly¬ 
mers only at the bottom wall. 

In order to produce a parabolic velocity profile inside 
the channel (e.g. along the z direction) a constant accel¬ 
eration A = Ax is applied to all fluid particles: 

mtfl luid = F{ 0 luid + m i A. (11) 

In turn, these fluid particle exchange moment with the 
chains, dragging them. Different values of A allow us to 
probe different dynamic regimes. 


2. BRUSH EQUILIBRIUM PROPERTIES 


Polymer brushes immersed in a liquid at equilibrium 
(i.e. in the absence of flow) have been thoroughly investi¬ 
gated in the past [55H4T] . Here we recall that the brush 
equilibrium conformation results from the balance be¬ 
tween configurational entropy, that tends to make chains 
visit the whole available space, and excluded volume in¬ 
teractions, avoiding contact between monomers. The 
brush can be properly described by its profile p(z), indi¬ 
cating the probability distributions of finding a monomer 
at distance z from the grafting wall. We show in Fig. [7] 
p(z) at fixed cr graft = 1-5 for several chain lengths 
(n = 20, 25, 30, 36,40,45). In line with previous theoreti¬ 
cal results [3S] and numerical studies [40132, the higher 
the volume fraction, the more step-like the density profile 
is. 

A definition of the brush height h in term of the 
first moment < z > of the density profile depends on 
the shape of p(z). Assuming that each chain is com¬ 
pletely elongated so that free ends are all at the max¬ 
imum distance from the grafting wall, the profile p(z) 
is step-like ( Alexander model) and the brush height can 
be defined as h = 2 < z >. A less rigid theory ( self- 
consistent field theory) hypothesizes that monomers of 
a same chain are distributed as a random walk, there¬ 
fore the free end position ranges uniformely over the 
whole brush thickness. In such a case the profile is 
parabolic, p(z) = p( 0 )(l — z 2 /h 2 ), and the brush height 

g 

h is h = — < z >. For the simulations discussed in the 

O 

manuscript we select n = 40, so that the evaluated p(z) is 
in between the step and the parabolic shapes (see Fig. [T]). 
In the following, we define h following the definition of a 
step distribution, e.g. 


0 /o°° zp(z)dz 

^ r pw 


2.1 Brush compression properties 


From simulations at equilibrium (A = 0) in which the 
brush is compressed by a flat surface kept at a fixed po¬ 
sition Sfi X we have estimated the Young’s modulus E as 
P 

E= ^ . In this expression, AS/S 0 = {S fix -S 0 )/S 0 

is the compression of the brush relative to the value So 
indicating the brush surface without any imposed con¬ 
finement and P is the pressure exerted by polymers on 
the pressing apparatus. In Fig. [ 8 ] we observe that E is in¬ 
creasing with the compression. We now notice that flow 
produces a compression of the brush (see Fig. [9]). Equat¬ 
ing the compression imposed at equilibrium with the one 
produced by the applied flow (shaded areas in Fig. [ 8 ]) we 











Conversion A <->■ Wi 



FIG. 7. Probability distributions of polymer brushes at 
equilibrium: fixed the grafting density, a gra ft = 1.5, we varied 
the degree of polymerization n. The case n = 40 corresponds 
to an almost step-like profile p(z). 


conclude that the brush elastic response is modified by a 
change in Wi. 



FIG. 8. The brush elastic response stored in the Young’s mod¬ 
ulus E. The modulus is defined as the coefficient between the 
relative brush compression and the applied pressure needed 
to produce it. The shaded areas indicate the ranges of brush 
surface compression corresponding to the three regions of flow 
inversion. 


3. THE WEISSENBERG NUMBER 

Instead of using A as a measure of the flow intensity we 
prefer to introduce an adimensional quantity, the Weis- 
senberg number Wi := tbrush/tfiow, where thrush should 
be a characteristic structural time of the unperturbed 
brush and tfi ow should provide an estimate of the typi¬ 
cal time scale associated to the flow. We evaluate tb rus h 
as the relaxation time obtained by the decay of the au¬ 
tocorrelation function, averaged over all distinct poly¬ 
mers, of the end-to-end vector amplitude R ee = ro — r n |. 


A 

Wi 

A 

Wi 

A 

Wi 

A 

Wi 

0.005 

45 

0.022 

208 

0.047 

409 

0.064 

612 

0.007 

64 

0.023 

218 

0.048 

413 

0.065 

616 

0.01 

91 

0.028 

275 

0.05 

418 

0.066 

627 

0.015 

137 

0.03 

297 

0.051 

420 

0.067 

642 

0.016 

147 

0.034 

340 

0.055 

430 

0.07 

737 

0.017 

158 

0.037 

375 

0.0565 

588 

0.075 

788 

0.018 

167 

0.04 

400 

0.575 

593 

0.08 

832 

0.019 

175 

0.043 

405 

0.058 

595 

0.1 

1044 

0.02 

187 

0.045 

406.7 

0.06 

599 



0.021 

200 

0.046 

407.4 

0.063 

607 




TABLE I. Tables of conversions between the constant acceler¬ 
ation A, imposed on each solvent particle to obtain a parabolic 
velocity profile, and the estimated a-dimensional Weissenberg 
number Wi. 



FIG. 9. The brush height h relative to the equilibrium value 
h eq as function of Wi. The brush height is defined in section 
2. For Wi < 100 the height keeps basically constant, while 
for Wi > 100 the flow has the effect to compress the brush. 


where r n and ro are the position vectors of, respec¬ 
tively, the last particle and the anchor of one chain in¬ 
side the brush. According to such definition, the relax¬ 
ation time tbrush assumes a dependency on the grafting 
density and on the degree of polymerization, thus repre¬ 
senting a time scale of the whole brush rather than an 
isolated single-chain property. The R ee auto-correlation 
function can be well represented by a stretched exponen¬ 
tial decay R ee ~ ; w here t 0 is the characteristic 

time scale and /3 the so-called stretching exponent. Thus 
tbrush = T[1 //3}t 0 /[3. For a gra ft = 1.5 and n = 40 we 
obtain tbrush = 1746. Regarding tfi ow , we choose the in¬ 
verse of the averaged shear rate 7 calculated as the ratio 
between the maximum flow velocity v ma x hi the channel 
and the corresponding z-coordinate R. Since our brush 
under flow experiences compression (see Fig. [ 9 ]) and poly¬ 
mers are set only at the bottom wall, both v max and R 
depend on the imposed acceleration A. 

In table Tab. |T] we list the correspondences between 
each value of A and its Wi. 
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4. VIDEOS 

We show in movies SI1 and SI2 the dynamics of a 
single chain inside the brush for two different Wi, one 
giving flow inversion (movie SI1) and the other with no 
flow inversion (movie SI2): the recursive motion is consti¬ 
tuted by an alternation of elongation-stretching-recoiling 
in both cases, thus it cannot be used as sufficient expla¬ 
nation for the backflow phenomenon. 

The motion of the whole brush under flow has been 
recorded in movies SI3 and SI4. In the first case no col¬ 
lective motion is detectable and no backward flow is ob¬ 
served in the velocity profile, while in the second case 
there is concurrence of a surface wave and of flow inver¬ 
sion. 


5. FINITE SIZE EFFECT 

We have investigated the influence of finite size effect 
increasing the channel length L x at fixed acceleration A: 
the wave persists. For a larger channel, 30 < L x < 40, 
the wave increases in wavelength (A = L x ) and decreases 
in frequency, as shown in Fig. [lOj such that the propaga¬ 
tion speed keeps constant. The wave amplitude decreases 
and, in accordance with the Taylor’s relation (see Letter, 
Eq. (1)), Vmi n decreases too. 

Another test has been attempted increasing, up till 
the double, both L x and L z (L x = 60 and L, = 100). 
For such a system we found again, for Wi = 464, a sur¬ 
face wave, with double amplitude and double wavelength 
(see movie SI5). However, investigating channels up to 
L x = 360, L y = 20 and L z = 150, we found that the 
amplitude of the surface wave does not keep on growing 
upon increasing L x : in fact, its value seems to reach a 
plateau. We interpret this non-linear behaviour as a con¬ 
sequence of the finite extensibility of the polymers, e.g. 
due to the higher and higher energy required to signifi¬ 
cantly compress and stretch the brush. 
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